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CARRIER PHASE ESTIMATION BASED ON SINGLE-AXIS CONSTANT MODULUS 
COST CRITERION AND BUSSGANG CRITERIA 



5 BACKGROUND OF THE INVENTION 

Field of the Invention 

The present invention relates to carrier phase and frequency detection in a receiver of a 
digital communications system. 



10 Cross-Reference to Related Applications 

This application claims the benefit of the filing date of U.S. provisional application no. 
60/248,941, filed on 1 1/15/00 as attorney docket no. 1043.002PROV. 



Description of the Related Art 

15 In many digital communications systems, a source generates digital information for 

transmission to multiple destination receivers. A transmitter processes the digital information 
into an encoded (e.g., error-correction encoded) and/or packetized stream of data. The stream of 
data is then divided into discrete blocks. Each of the blocks is mapped onto a corresponding one 
of a sequence of code or symbol values ("symbols") chosen from a pre-defined alphabet A, and 

20 generated with a period T s , sometimes referred to as the "baud" period. Symbols may be 
modulated by an analog, e.g., radio frequency (RF), carrier, in amplitude, phase, and/or 
frequency prior to physical transmission through the communication medium. Many methods of 
mapping exist and are well known in the art, and these pre-defined alphabets are generated based 
on certain criteria. For example, data may be mapped into symbols of a complex data stream as 

25 pairs of in-phase (I) and quadrature phase (Q) component values that are subsequently 
modulated with an RF carrier. 

A receiver performs several functions to demodulate and decode a received signal. 
Receiver functions include, for example, tuning and RF demodulation of the received signal to 
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an intermediate frequency (IF) signal, synchronization with the RF carrier, equalization, symbol 
detection, and decoding. 

FIG. 1 shows a typical prior art communication system 100 that may be employed for 
transmission and reception of digital television signals. Communication system 100 comprises 
5 transmitter 101 transferring signals through transmission medium 102 to receiver 103. 

Transmitter 101 comprises digital encoding system 111, premodulator/pulse shaper 1 12, radio 
frequency (RF) upconverter 113, carrier oscillator 115, and transmit antenna 1 14. Transmitter 
101 receives user data from information source 1 10 (such as video, audio, and/or computer files) 
coupled to digital encoding system 111. Digital encoding system 1 1 1 may provide analog-to- 

10 digital (A/D) conversion, error-correction encoding, and/or bit-to-symbol mapping to generate a 
sequence of symbols selected from a predetermined alphabet. For example, data may be mapped 
into a complex-valued signal stream with pairs of in-phase (I) and quadrature phase (Q) 
components. Digital encoding system 1 1 1 provides the symbols to pre-modulator/pulse shaper 
112. Pre-modulator/pulse shaper 1 12 modifies the symbols for the particular type of modulation, 

15 and may include a filter for pulse shaping of the symbols. The signal generated by pulse shaper 
1 12 is provided to RF upconverter 113 which i) uses the signal to modulate a complex radio 
frequency (RF) carrier provided by carrier oscillator 115, and ii) amplifies and filters the signal. 
The modulated and amplified RF carrier is then emitted into the transmission medium 102 as an 
RF signal via transmit antenna 1 14. 

20 Various modulation techniques, such as quadrature amplitude modulation (QAM), phase- 

shift keyed (PSK) modulation, or vestigial sideband (VSB) modulation are known in the art of 
communications to modulate the carrier. For example, modulation formats such as VSB are 
common formats used for transmission of digital television signals in accordance with, for 
example, the ATSC standard for digital television, "ATSC Digital Television Standard," Doc. 

25 A/53, Sept. 1995. 

For these modulation techniques, a quadrature oscillator may be employed with a 
complex RF upconverter to form a modulator. The I signal component modulates the cosine 
component generated by the oscillator and the Q signal component modulates the sine 
component of the oscillator. VSB modulation is a form of single-sideband modulation in which 
30 the redundant sideband of a real- valued signal is removed in full by filtering, except for a small 
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vestige of the sideband. For complex VSB modulation, a complex signal is formed with the Q 
component being the Hilbert transform of the I component (however, the Q-component thus 
contains no additional user information). 

FIG. 2 shows complex VSB modulation implemented by premodulator/pulse shaper 1 12, 
5 carrier oscillator 115, and RF upconverter 113 of transmitter 101 inFIG. 1. The signal from 
digital encoding system 1 1 1 is split into two paths: an upper path and a lower path. The signal in 
the upper path is filtered with Hilbert filter 201 and shifted in phase by 90° with multiplier 204, 
creating the Q (imaginary) component of the VSB signal. The signal in the lower path is applied 
to delay 202, creating a delayed version of the lower signal corresponding to the I (real) 
10 component of the VSB signal. The signals from the upper and lower paths are combined in 
adder 205 to form the complex VSB signal. The complex VSB signal excites pulse shaper 203 
to generate pulses having a shape optimized for detection at the receiver. Up-converter 1 13 then 
up-converts the output signal of pulse shaper 203 to the RF carrier frequency with a complex RF 
modulator using the complex carrier generated by oscillator 115. 

1 5 The modulated carrier signal transmitted through the medium 102 (which may be, e.g., 

wire, optical fiber, atmosphere, space, etc.) comprises a series of analog pulses, each analog 
pulse being amplitude and/or phase modulated by a corresponding symbol in the sequence. The 
pulse shape used typically extends many symbol periods in time. This introduces the possibility 
of adjacent pulses corrupting each other, a phenomenon known as inter-symbol interference 

20 (ISI). 

Returning to FIG. 1, receiver 103 includes antenna 120 receiving the signal from the 
medium 102, complex demodulator and sampler 121, timing recovery module 122, equalizer 
123, and carrier recovery module 124. Carrier recovery module 124 includes reference 
generator 126 and phase detector 125. Complex demodulator and sampler 121 translates the 

25 received signal from RF to intermediate frequency (IF), and performs complex demodulation of 
the received signal at IF to near passband employing the locally generated reference for the 
cairier signal. Complex demodulator and sampler 121 also samples the signal based on an 
estimate of the symbol period. Timing recovery module 122 estimates the symbol timing period 
T s , and this estimate may be fed back to complex demodulator and sampler 121 to adjust the 

30 sampling rate (e.g., via a sampling phase error). Timing recovery thus synchronizes sampling 
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instances to the top-dead-center of the pulse shapes, and then tracks variations in the detected 
period and phase of T s . 

Equalizer 123 applies equalization to the received samples, e.g. to correct for ISI, and 
may generate a cost error term used by timing recovery module 122 to adjust its estimate of the 
5 symbol timing period T s . Carrier recovery module 124 generates estimates for the difference in 
frequency and phase (collectively referred to as angle 0) of the carrier used to modulate the 
symbols and the locally generated reference used for demodulation. From the estimate of angle 
0, carrier recovery module 124 adjusts the reference for complex demodulation to adjust the 
received signal to precise baseband. A detector 150, typically implemented with a slicer, 
10 examines each sample to generate either a soft or hard decision for the symbol that corresponds 
to the sample. After symbol detection, a decoder 151 reconstructs the transmitted data from the 
symbol sequence. 

Many methods exist in the art for achieving synchronization of sample timing and/or 
carrier recovery. One method uses a separate pilot tone or preamble pattern in phase with the 

15 modulation process that is transmitted in addition to the information-bearing signal. The 

receiver derives synchronization measures from the separately transmitted information e.g., from 
the pilot tone. However, including a separate signal for synchronization reduces channel 
throughput (and uses both extra power and bandwidth) for pilot tone transmission and reception. 
Consequently, many applications use blind recovery techniques in which information for 

20 synchronization is derived from the received signal itself without the aid of side information, 
such as a pilot tone. 

A receiver generally requires accurate knowledge of the received carrier signal's 
frequency and the carrier signal's phase offset (i.e., angle 0). Estimates of carrier frequency and 
phase offset are required to properly recover I and Q components of a passband or near-baseband 
25 signal and adjust the received signal to precise baseband for symbol detection. The estimation of 
carrier frequency and phase may be generated using samples obtained prior to equalization by 
equalizer 123, or equalized samples from equalizer 123. 

A two-step procedure is often used in the prior art to accurately estimate the carrier 
frequency and phase offset at a given moment. The first step derives a coarse estimate of the 
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frequency and phase offset, which reference generator 126 may generate directly from the 
samples. The second step derives a fine (i.e., an adaptive) estimate, and a phase tracking loop is 
generally employed including reference generator 126 and phase detector 125. The fine 
estimate tracks the phase jitter introduced into the received signal by 1) time- varying channel 
5 impairments and 2) phase noise introduced by low-cost tuner oscillators. Carrier estimation may 
be performed in discrete time, in continuous time, or in a combination of both discrete and 
continuous time. 

The carrier tracking loop may use an estimated error as a control signal to the reference 
generator to minimize rotation between I and Q components of a locally generated estimate of 
10 the carrier signal ("reference") and the received carrier signal. Alternatively, the estimated error 
may be processed to generate the error angle G that is then used to de-rotate the demodulated 
signal Phase detector 125 generates the measure (i.e., the error angle 9) between the reference 
q and received carrier signal. The reference may be generated from a separate synchronization 
signal at the transmitter (termed "trained estimation") or derived from the information signal 
fyl 5 itself (termed "blind estimation"). Trained estimation may be based on a continuous wave, such 
q as a pure tone signal or a pseudo-random (PR) sequence of digital symbols, generated by the 
jzf transmitter. For systems using blind estimation, the reference is not known to the receiver, and 
s two typical techniques employed in the prior art are "power of N" and "decision directed" carrier 
m recovery. 

\20 In power-of-N carrier recovery, the received signal is raised to the power N to create a 

r: strong spectral line at frequency 27rNf c , where N is positive and f c is the carrier frequency. The 
phase detector compares the phase of the signal raised to the power N to an initial estimated 
phase. The output (error) signal of the phase detector is loop filtered and applied to a voltage- 
controlled oscillator (VCO) with natural frequency 27rNf c . The output of the VCO corresponds to 
25 the estimated phase of the transmitted signal. The estimated offset is fed back to the phase 

detector and compared to the signal raised to the power N. The procedure is repeated until the 
error signal at the output of the phase detector is minimized. Power-of-N carrier recovery 
generally produces a noisy estimate of the carrier phase and frequency. 

Decision-directed carrier recovery employs a phase detector that measures the angle 6 
30 between the received and (estimated) reference I and Q components that correspond to a nearest 

^ Attorney Docket No. : 1 043 , 002 



predefined symbol or "alphabet" member. Since each transmitted character is formed from 
known I and Q component values in amplitude and phase, the reference I and Q components are 
derived from the estimate (receiver's guess) of which symbol is received and under study. The 
output of the phase detector is coupled to a loop filter that controls the VCO generating the 
reference I and Q components. The output of the VCO is fed back to the phase detector to 
complete the phase tracking loop. Decision-directed carrier recovery is typically employed 
when the overall system Signal-to-Noise Ratio (SNR) is high. 

Trained estimation is similar to blind estimation, except that the receiver uses the 
received reference to reduce error from estimation of the reference. For example, in carrier 
estimation based on a pilot tone, the estimate is derived in a manner similar to that of the power- 
of-N blind carrier estimate, but a pure tone signal of arbitrary frequency sent by the transmitter 
and known at the receiver aids in carrier estimation. For the ATSC standard for broadcast of 
digital television in the U.S., the transmitter embeds a single pilot tone at the band edge of the 
data spectrum. Mean square error estimation is comparable to the decision-directed blind 
estimation, but the angle between the components is computed by minimizing squared Euclidean 
distance between I and Q components of an estimated symbol sequence and a sequence 
generated by the transmitter (i.e., the training reference signal). 

A receiver also generally applies equalization to the sample sequence prior to forming 
hard decisions for symbols. Equalization is used to suppress the effects of ISI, caused by 
phenomena such as i) residual timing error, ii) pulse shape/multipath distortion from the 
propagation channel, and/or iii) approximation of the ideal transmit and receive filters for ease of 
implementation. Adaptive equalizers may also use blind recovery techniques to derive tap 
coefficients for the equalizer filters. 

One such blind cost criterion employed for equalization is the constant modulus (CM) 
criterion. The stochastic gradient descent of the CM criterion is known as the Constant Modulus 
Algorithm (CMA). The CMA algorithm is described in an article by D.N. Godard entitled "Self- 
Recovering Equalization in Two-Dimensional Data Communication Systems," IEEE 
Transactions on Communications, vol. 28, no. 11, pp. 1867-1875, Oct. 1980, which is 
incorporated herein by reference. The CM criterion and CMA algorithm were further developed 
to de-couple equalization and carrier recovery functions in a receiver. Such use of the CM 
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criterion and CMA algorithm for equalization is described in J.R. Treichler et al., "A New 
Approach to Multipath Correction of Constant Modulus Signals," IEEE Transactions on 
Acoustics, Speech, and Signal Processing , vol. ASSP-31, no. 2, April 1993, which is 
incorporated herein by reference. Systems that use such CMA algorithm for adaptive 
5 equalization, such as that described in U.S. Patent No. 5,872,815 to Strolle et al., do not employ 
the CM criterion for timing or carrier phase recovery. 

The CM criterion penalizes the deviation of the dispersion of the magnitude squared of 
the received signal from a pre-calculated constant referred to as the "dispersion constant" or the 
"Godard radius." FIGs. 3A and 3B illustrate that the constant modulus criterion is based on 
10 deviation from a "radius" about the origin of, for example, a source constellation. FIG. 3 A 
shows a radius 301 of an 8-PSK (phase-shift keyed) constellation plotted for real (e.g., Re or I) 
and imaginary (e.g., Im or Q) components. In FIG. 3A, each point (symbol) lies on the circle 
^ defined by this radius (termed a constant modulus system), and the CM criterion penalizes 
tff distance of a received sample (e.g., sample 302) from this circle. Even though the constellation 
f||l 5 may rotate, the constellation remains on the circle, and so applying a CM criterion to this 
^ constellation does not penalize spatial rotation of the constellation due to residual carrier offset. 
Q FIG. 3B shows a 16-QAM (quadrature amplitude modulation) constellation plotted for real and 
s imaginary components. In FIG. 3B, groups of points (symbols) lie on corresponding concentric 
% circles 311,312, and 313. The CM criterion determines a radius 3 1 4 of circle 315, which is not 
CEO necessarily a radius of one of the concentric circles 31 1, 312, and 313 (termed non-constant 
□ modulus), as a "common" radial distance from the origin for the points of the constellation. As 
s "~ with the constellation of FIG. 3 A, the CM criterion penalizes distance of a received sample (e.g., 
sample 303) from this circle 315. 

The CM criterion defines a cost function J C m that may be expressed as given in equation 

25 (1): 

JcM = E[(p 2 -|y n (x,g)| 2 ) 2 ] (1) 

where E[«] denotes the expected value, p 2 is the dispersion constant (Godard radius), y n (x,g) is 
the discrete value that represents the sampled signal, x represents the timing (sampling) phase, 
and g represents the equalizer taps introduced to suppress the ISI. The dispersion constant p 2 is 
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a quantity that can be determined from the type of modulation employed (e.g., QAM, BPSK, 
etc.). The dispersion constant p 2 may be derived by calculation, by experiment, or by a 
combination of both for a particular implementation. For real-valued source, such as VSB, the 
CM criterion, and its stochastic gradient, may be modified by taking the real part of y n (x,g) in 
5 equation (1). The modified CM criterion is referred to as the single-axis (SA) CM criterion, and 
is given in equation (2). 

Jsa-cm = E[(p 2 - Re{y n (T,g, 9)} 2 ) 2 ] (2) 

where Re{«} denotes the real-part extraction. 

Given a defined cost function, the gradient of the cost function may be derived. The 
10 stochastic gradient is an approximation of the true gradient that is calculated by taking the 
derivative of the cost function without taking the expected value. For example, the stochastic 
gradient of the single-axis (S A) CM criterion is known as the SA-CMA and is derived by taking 
,Q the derivative of equation (2) with respect to the variable of interest. Once the derivative is 
2j1 calculated, an error term may be defined that tends to drive the stochastic gradient to a relative 
[Li 5 minimum. For timing recovery, derivation of such stochastic gradient is described in U.S. Patent 
p application no. 09/761,303, entitled "Blind Cost Criterion Timing Recovery," filed on Jan. 17, 
* : 200 1 , by the inventors herein and is incorporated herein by reference. 

m SUMMARY OF THE INVENTION 

u\ The present invention relates to carrier recovery using carrier frequency phase error 

0 estimation employing a single-axis constant modulus (SA-CM) criterion to generate the carrier 
frequency phase error. The carrier phase error is used to adjust a received and demodulated 
signal from near to precise baseband. For one implementation, the estimate is derived in 
accordance with the stochastic gradient of the SA-CM criterion via the SA-CM algorithm (SA- 
CMA). The carrier phase error is then used to adjust the carrier frequency and phase of the 
25 received and demodulated signal toward the phase of the carrier used to modulate the transmitted 
symbols, driving the carrier phase error to a constant, such as zero. 

In accordance with an exemplary embodiment of the present invention, carrier recovery 
in a demodulated signal formed from a data-modulated carrier i) generates an estimate of an 
angle between the carrier and the locally generated reference from the signal and based on a 
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stochastic gradient of a single-axis cost function, the cost function being selected from a set of 
Bussgang-class cost functions; and ii) adjusts at least one of the frequency and phase of the 
demodulated signal based on the angle such that the magnitude of the angle is driven to a 
predetermined value. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Other aspects, features, and advantages of the present invention will become more fully 
apparent from the following detailed description, the appended claims, and the accompanying 
drawings in which: 

FIG. 1 shows a typical prior art communication system that may be employed for 
transmission of digital television signals; 

FIG. 2 shows a prior art vestigial sideband (VSB) modulator that may be employed by the 
transmitter of FIG. 1; 

FIG. 3 A illustrates a constant modulus criterion based on determining a radius about the 
origin of an 8-PSK source constellation; 

FIG. 3B illustrates a constant modulus criterion based on determining a radius about the 
origin of a 16-QAM source constellation; 

FIG. 4 shows a receiver of a digital communication system that employs carrier recovery 
in accordance with embodiments of the present invention based on a single-axis, blind cost 
criterion; 

FIG. 5 shows a block diagram of the functional components of the detector of FIG. 4 
including forward filter and feedback equalizer/carrier loop (FECL) module; 

FIG. 6 shows an exemplary implementation of the FECL module of FIG. 5 in which the 
feedback filter operates with signals at the baseband; 

FIG. 7 shows an exemplary implementation of the FECL module of FIG. 5 in which the 
feedback filter operates with signals at the passband; 

FIG. 8 shows an exemplary embodiment of the carrier tracking loop of FIG. 5; 
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FIG. 9 shows an exemplary implementation for a portion of the carrier tracking loop 
operating in accordance with a first configuration of the FECL module of FIG. 6; 

FIG. 10 shows an exemplary implementation for a portion of the carrier tracking loop 
when the carrier phase error term is calculated externally to the carrier tracking loop; and 

FIG. 1 1 shows an exemplary carrier tracking loop employing adaptive switching between 
a S A-CM phase detector and other types of phase detectors based on relative quality of the 
received signal. 

DETAILED DESCRIPTION 

FIG. 4 shows a receiver section 400 of a digital communication system that employs 
residual carrier phase recovery and synchronization (termed herein as carrier recovery) based on 
the single-axis blind cost criterion in accordance with embodiments of the present invention. 
Carrier recovery determines an estimate of the angle 0 representing the phase error between the 
transmit carrier signal and the locally generated carrier signal, and uses the estimate to minimize 
the phase error, such as by driving the phase error to zero. The preferred embodiments described 
herein employ the single-axis constant modulus (SA-CM) cost function and a carrier phase error 
term derived from its stochastic gradient (the stochastic gradient descent of the SA-CM cost 
function) to generate an estimate of the angle 9. Carrier recovery in accordance with 
embodiments of the present invention generates a carrier phase error employing the SA-CM 
error term and a derivative of the demodulated signal with respect to the angle 0. Embodiments 
of the present invention tend to adjust the residual carrier frequency/phase of the received signal 
in a manner opposite to that of the error term, and thus drive the carrier phase error, and the 
angle 0, to zero. 

Receiver 400 receives modulated signals from transmission medium 402, and receiver 
400 may be employed in a digital communication system broadcasting digital television signals 
in accordance with, for example, the ATSC standard for digital television, "ATSC Digital 
Television Standard," Doc. A/53, Sept. 1995. Consequently, the preferred embodiments of the 
present invention operate on received signals having complex VSB modulation format. VSB is 
the modulation format adopted for terrestrial broadcast of digital television signals in the United 
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States by the ATSC standard. 

Figure 4 shows the signal path of I and Q components (and processing) as a single line, 
though they may be conveyed and/or processed through separate paths. In general, signals 
within receiver 400 may be complex- valued, though for some applications only I components are 
5 used. In addition, while the figures illustrate carrier recovery in a digital television receiver, 
many components in addition to those shown may exist, but these additional components are not 
necessary to understanding the present invention and are omitted for simplicity. 

Returning to FIG. 4, receiver 400 comprises receive antenna 420, tuner 421, 
demodulator/downconverter 422, timing recovery module 423, detector 424, and decoder 425. 
10 Antenna 420 receives the RF signal from the transmission medium 402. The received RF signal 
may have added phase and amplitude noise, doppler shift, multipath effects, distortion, 
dispersion, and/or other effects typically associated with the particular transmission medium. 

O Tuner 421 isolates the RF signal of interest received through antenna 420 from medium 

if) 402 and translates the signal, in frequency, to an intermediary frequency (IF) signal ("passband 
LS 5 signal"). Many different techniques exist in the art for demodulating a modulated RF carrier to 
□ provide a near-baseband or passband signal, and some of these systems may not necessarily 
m employ digital signal processing techniques. One skilled in the art would realize that the process 
;L. of receiving, isolating, and demodulating the transmitted signal is described herein as an aide to 
61 the description of carrier recovery in accordance with exemplary embodiments of the present 
30 invention. 

M= Demodulator/downconverter 422 processes the passband signal as follows. The received 

signal is sampled via an analog-to-digital (A/D) converter 450, creating a digital sample 
sequence. A/D converter 450 may generate sample values at a rate greater than the symbol 
period. As would be apparent to one skilled in the art, decimation and/or interpolation may be 

25 employed at various points within receiver section 403 to generate sample sequences at different 
sampling rates and different sampling phases. Demodulator/downconverter 422 may apply pre- 
filter/power detector 45 1 to the digital sample sequence to reduce noise effects and to generate a 
power level estimate for automatic gain control (AGC). The gain of the near-baseband signal is 
adjusted via an automatic gain control (AGC) circuit 452. The gain of the IF signal and/or RF 

30 signal may also be controlled via AGC circuits implemented similarly to that of AGC circuit 
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452. AGC circuit 452 controls the dynamic range of the near-baseband signal to within a 
predetermined range. 

Complex demodulator 453 demodulates the received signal to a passband signal or to a 
near-baseband signal, and splits the demodulated signal into a complex sample sequence having 
5 in-phase (I) and quadrature phase (Q) components. Complex demodulator 453 may employ an 
oscillator generating sine and cosine components at a frequency and phase based on information 
derived by phase-locked loop/carrier recovery (PLL/CR) module 462 of detector 424 operating 
in accordance with an exemplary embodiment of the present invention, as described 
subsequently. The near-baseband signal is not at precise baseband due to sampling clock 
10 inaccuracies, oscillator phase noise, or doppler shift. Timing recovery module 423 aligns or 
otherwise synchronizes the period T and phase t for sampling of the received signal to the 
symbol period (baud rate) T s and phase x s of I and Q signal components of the transmitted signal, 
Q respectively. 

*0 The sequence of samples is then provided to detector 424, which generates an estimate 

H15 for each sample corresponding to a transmitted symbol. Detector 424 filters the signal with 
y matched filter 460 whose frequency response characteristic matches the frequency response of 
Qri the transmitter's pulse shaper. For example, if the pulse shaper applies a filter corresponding to 
q a square root-raised cosine function, the matched filter also applies a square root-raised cosine to 
~f i) create a Nyquist pulse at the matched filter output and ii) maximize signal-to-noise ratio 
lo (SNR). 

N 5 Since transmission medium 402 may distort the signal waveform passing through it, each 

sample from timing recovery module 423 may include ISI effects of channel multipath. 
Equalizer 461 of detector 424 applies equalization to the sequence of samples from timing 
recovery module 423. Equalizer 461 may be implemented with an adaptive, multi-tap filter in a 

25 linear or non-linear architecture and with a finite or infinite impulse response. Equalizer 461 

adjusts the sample values to remove distortion effects and restore the "zero ISI" characteristics of 
the sample sequence. For highly impaired media, such as media experiencing long echo 
multipath, decision feedback equalization (DFE) may be preferred over a linear filter. A decision 
feedback equalizer acts to additively cancel ISI. 
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PLL/CR 462 provides residual carrier synchronization based on a single-axis constant 
modulus (SA-CM) cost criterion in accordance with embodiments of the present invention. 
Carrier recovery translates the (sampled) signal from near to precise baseband by compensating 
for a residual offset between the carrier frequency of the transmitted signal and frequency for 
demodulation selected by tuner 421 and/or complex demodulator 453. Symbol decision module 
463 generates a hard decision for each symbol from the corresponding equalized/compensated 
samples. While shown as part of detector 424, either detector 424 or decoder 425 may employ a 
symbol decision module 463, which may be typically implemented with a slicer, nearest element 
quantizer, maximum likelihood (ML) detector, or maximum a posteriori (MAP) detector. 

Decoder 425 processes the sequence of symbol decisions to reverse the various levels of 
encoding (symbol-to-data mapping, error correction, re-formatting, and/or depacketizing) to 
reconstruct the original user information. Re-formatting may include decoding in accordance 
with a high-definition digital television standard (e.g., MPEG-1, MPEG-2). 

In accordance with embodiments of the present invention, carrier recovery and 
synchronization by, for example, PLL/CR 462 operates "blindly". Thus carrier recovery and 
synchronization i) operate without using a reference signal, ii) includes a phase detector 
implementing adaptive carrier phase frequency estimation, and iii) the carrier loop generates the 
precise carrier frequency. The preferred embodiments of the present invention are described 
herein with respect to adaptive carrier-phase estimation employing the single-axis, constant 
modulus (SA-CM) criterion and its stochastic gradient descent (SA-CMA) providing the phase 
error terms csa-cm- 

The SA-CM criterion may be defined for VSB signals, such as described in a paper by 
Shah et al, "Global convergence of a single-axis constant modulus algorithm," Proceedings of 
the Statistical Signal and Array Processing Workshon. Pocono Manor, PA, August 2000, which 
is incorporated herein by reference. SA-CM for VSB signals and blind cost error terms 
generated from the stochastic gradient of SA-CM criterion are also described in U.S. Patent 
application no. 09/100,705, entitled "A REDUCED COMPLEXITY BLIND EQUALIZER FOR 
MULTI-MODE SIGNALING," filed on June 19, 1998. The SA-CM criterion is typically based 
on real-part extraction of the output signal after equalization, and may be used for both real- or 
complex-valued equalizer coefficients operating on either baseband or passband signals. 
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In accordance with embodiments of the present invention, a carrier tracking loop of 
PLL/CR 462 includes SA-CM based phase detection. For this case, the signal considered is a 
discrete signal present in detector 424 that is a function of the carrier phase angle 9. Notice that 
contrary to prior art equation (1), equation (2) is a function of the phase of the sampled signal y n 
The cost function Jsa-cm can therefore be jointly optimized with respect any of the variables x,g, 
or 9. Defining sample y n (9) equivalently as y n (x,g,9) to indicate that the sample (a complex- 
valued signal) is a function of phase angle 9 which is the parameter of interest, then the SA-CM 
criterion is given by equation (3): 

JcM = E[(p 2 -Re{y n (9)} 2 ) 2 ]. (3) 

where p 2 is the Godard radius, Re{#} denotes real-part extraction, and E<» denotes the expected 
value. The single-axis gradient descent may be determined by taking the expected value of the 
derivative of the cost function with respect to y„(8). The inventors have determined that an 
approximation of the gradient may be calculated without taking the expectation operation. Thus, 
the single-axis stochastic gradient at discrete time [n], is given by equation (4): 

dJ CM /d9 = 4 e SA -cM[n] dy n (9)/d9 (4) 

where the SA-CM phase error term at discrete time [n], e S A-CM[n], is given by equation (5): 

e S A-cM[n]- Re{y n (9)}(p 2 - Re{y n (9)} 2 ). (5) 

The value p 2 is a known constant for the type of modulation. The constant p 2 is calculated as 
the ratio of the fourth to second central source moments. For example, p 2 is 37 when VSB-8 
modulation format is used with symbol values +/- 1, +/-3, .. ., +/-7. 

The derivative of y n (9) with respect to the angle 9 may also be approximated as DT[n] for 
the purpose of calculating quantities from the sample sequence. The stochastic gradient SA- 
CMA employed by the phase detector is then expressed as in equation (6): 

dJ CM /d9 = 4 e SA -CM[n] DT[n], (6) 

where the SA-CM phase error term is the same as equation (5), and where DT[n]=dy n (9)/d8 is 
the exact or approximate derivative of y n (9) with respect to the angle 9. SA-CM phase detection 
computes the product of the two quantities e SA -CM[n] and DT[n]. The "n" shown as a subscript 
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and shown in brackets with respect to a variable is equivalent and indicates a digital sample 
value at discrete time n. It may be shown that operating in accordance with equation (7) is 
operation as a phase detector. 

To track phase angle 9[n], a second order phase lock loop with SA-CM as the phase 
detector may use the update rule as given in equation (7): 

9 [n] = 6 [n-1] + 4 F(z) e SA -cM[n] DT[n] (7) 

where 6 [n] is defined as an arbitrary rotation in phase (phase error) and where F(z) is a digital 
filter of the form given in equation (8): 

F(z)=7(Yi +Y2/(l-z ( - 1) )), (8) 

where y, yi and y 2 are programmable gains selected for the particular implementation. The first 
order implementation is referred to the case where j2 =0. 

FIG. 5 is a block diagram showing the functional components of detector 424 including 
forward filter 501 and feedback equalizer/carrier loop (FECL) module 502 operating in 
accordance with the present invention. FECL module 502 includes carrier tracking loop 503, 
decision device 504, and feedback filter 505. Carrier tracking loop 503 generates estimates of 
the carrier phase error, denoted as the angle 8, in accordance with embodiments of the present 
invention. 

Decision device 504 is included in symbol decision module 463 (FIG. 4) and is employed 
by detector 424 to create a "hard" decision for an alphabet element (symbol) corresponding to an 
input sample. Decision device 504 may be typically implemented with a slicer. Several methods 
may be employed by a slicer to generate a hard decision for systems using VSB modulation. 
These methods include: (i) calculating the I component of the sample to the nearest (most likely) 
alphabet element and not processing the Q component of the sample; (ii) quantizing the I 
component of the sample and quantizing the Q component so that the phase angle of the complex 
input sample is unchanged; and (iii) quantizing the I component of the sample and generating the 
Q component of the sample using a Hilbert transform of the I component of the sample. 

Alternatively equalizer 461 of FIG. 4 may be an adaptive, multi-tap feedback equalizer 
using finite impulse response forward and feedback filters to realize an infinite impulse response 
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(IIR) equalizing filter. Forward filter 501, decision device 504, and feedback filter 505 are 
configured such that the output samples of the equalizer are filtered by feedback filter 505 and 
added to the output of forward filter 501. Equalizer 461 may be implemented with forward filter 
501 and feedback filter 505 if equalizer 461 is a linear IIR equalizer (i.e., no hard or soft decision 
5 of decision device 504 is input to feedback filter 505). Alternatively, equalizer 461 may be 
implemented with forward filter 501, feedback filter 505, and decision device 504 if equalizer 
461 is a feedback equalizer (i.e., a hard or soft decision of decision device 504 is input to 
feedback filter 505). A preferred embodiment of the present invention allows for switching 
between linear and decision feedback equalization operation. 

10 Methods to generate tap values for the equalizer 461 are known in the art. For example, 

the CM criterion and CMA algorithm for equalization are described in J.R. Treichler et al., "A 
New Approach to Multipath Correction of Constant Modulus Signals," IEEE Transactions on 
□ Acoustics, Speech, and Signal Processing , vol. ASSP-3 1, no. 2, April 1993; and in U.S. Patent 
^ No. 5,872,815 to Strolle et al. As referenced above, U.S. Patent application no. 09/100,705 
fUS describes use of the SA-CM criterion and SA-CMA for adaptive equalization. Some 
p embodiments of the present invention may employ the techniques of U.S. Patent application no. 
y 09/100,705 to generate tap values for forward filter 501 and feedback filter 505 when employed 
B as a DFE. 

£9 Returning to FIG. 5, interconnection of carrier tracking loop 503, decision device 504, 

gp and feedback filter 505 of FIG. 5 depends on when "de-rotation" of the signal samples (i.e., 
} := -; translation of samples to precise baseband) occurs in the implementation. One skilled in the art 
will realize that each implementation is a consequence of the model employed to describe 
perturbations on the channel phase. FIGs. 6 and 7 show implementations 502a and 502b, 
respectively, of FECL module 502 shown in FIG. 5 that receive the output samples of forward 
25 filter 501. FIG. 6 shows an implementation of FECL module 502a in which feedback filter 505 
operates with signals at the baseband. FIG. 7 shows an implementation of FECL module 502b in 
which feedback filter 505 operates with signals at the passband. 

Referring to FIG. 6, the complex (FIR) signal (sample stream from forward filter 501) is 
applied to phase rotator 601. Phase rotator 601 operates on the complex (I and Q) components, 
30 denoted y n (0Mn(6) + jQ n (Q), and translates the samples to baseband by computing the product 
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y n (6)e" J [n] , where 0[n] is an estimate at discrete time n of the carrier phase error (offset angle 0). 
The rotated signal is added combined with the output of the feedback filter 505 by real combiner 
602. The output signal of combiner 602 is applied to decision device 504, which maps the 
output signal to the nearest alphabet symbol. When switch 603 is in position #1, the input to 
feedback filter 505 is the output of combiner 602 (e.g., when implementing a linear equalizer). 
When switch 603 is in position #2, the input to feedback filter 505 is the output of decision 
device 504 (e.g., when implementing a DFE). 

Carrier tracking loop 503 receives either the output decisions of decision device 504 
(e.g., the output symbol of the slicer), the output signal from combiner 602, or the output of the 
phase rotator 601, depending on which mode of operation. Modes of operation include decision- 
directed (DD) and blind (using SA-CMA phase detection as described herein) modes. For DD 
mode carrier phase estimation, carrier tracking loop 503 receives both the signal from combiner 
602 and output decisions of decision device 504, and from these inputs generates the estimate for 
angle 0. For blind mode carrier phase estimation, carrier tracking loop 503 includes a blind 
phase detector to calculate the value for carrier phase error given in equation (7), which value is 
then (approximately) integrated to give the output angle 0. 

Embodiments of the present invention operating in the baseband such as shown in FIG. 6 
may generate an estimate of 0[n] using one of at least three baseband configurations. In a first 
exemplary configuration of FECL module 502a, carrier tracking loop 503 generates an estimate 
of 0[n] based on the input Re{y n (0)e" j9[n] } provided by the I component output from phase rotator 
601. In a second exemplary configuration of FECL module 502, carrier tracking loop 503 
generates an estimate of 0[n] based on the input Re{y n (0)e" je[n] - w[n]}, where w[n] is the output 
of the feedback filter. In this second exemplary configuration, carrier phase correction is based 
on the 1 component output from combiner 602. In a third exemplary configuration of FECL 
module 502, carrier tracking loop 503 generates an estimate of 0[n] based on the output of 
decision device 504 and combiner. 

As described previously, SA-CM phase detection computes carrier phase error with the 
stochastic gradient dJcivi/dO as given in equation (6) and repeated for convenience as follows: 

dJ CM /d0 = 4e SA -cM[n] DT[n], (6) 
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with e SA -CM[n]= (Re{y n (9)} 2 - p 2 )Re{y n (9)} and DT[n]=d(y n (9))/d9 are as defined previously. 

At discrete time n, the SA-CM phase detector of PLL/CR 462 operates on the previous 
estimate of 0[n-l] and on the inputs Re{y n (9)} and Im{y n (9)}. The signal Re{y n (9)} is used to 
compute e S A-cM[n] by squaring Re{y n (9)}, subtracting p 2 and multiplying the difference by 
5 Re{y n (0)}. The phase detector signal is formed by multiplying e S A-cM[n] with DT[n] to drive a 
loop filter that adjusts the carrier phase of the sample sequence. 

However, the data signal input to the carrier loop may be generated at different locations 
of the configurations of FIG. 6, and so values for e S A-CM[n] and DT[n] vary depending on the 
particular configuration. The location where the signal for the phase detector is generate depends 
10 implicitly of the model used for the channel. 

FIG. 8 shows an exemplary embodiment of carrier tracking loop 503 with SA-CM phase 
p detection. Signal data provided to carrier tracking loop 503 (e.g., generated from various paths 

as described with respect to FIG. 6 and 7) is applied to SA-CM phase detector 801 . SA-CM 
fu phase detector 801 generates a carrier phase error, which is an estimate of the derivative of the 
Hi 5 input signal with respect to the measure (in this case, phase angle 9[n]) scaled by the blind cost 
~ criterion (i.e., SA-CMA) error term. The carrier phase error is filtered by (discrete-time) loop 
^ filter 802 driven by the parameters yi and y 2 , and scaled by a constant y in scaler gain 803. To 
£B generate an estimate of phase angle 9[n], the output of scaler gain 803 is integrated. Combiner 
j™j 804 and delay 805 implement an approximation to the integral to generate the estimate of phase 
jlo angle 9[n]. Module 806 computes <j m by generating the cosine and sine of phase angle 9[n]. 

For the first configuration of FECL module 502a shown in FIG. 6, the signal data is the 
output signal from phase rotator 601 . The phase error term e SA -cM[n] of SA-CM detector 801 is 
computed with the data signal y n (9)e" je[n " I] using the previously generated value of phase estimate 
9[n-l]. The contribution DT[n] is expressed as in equation (9): 

25 DT[n] = Re{y n (9)e- je[n] }= Re{y n (9)}cos(9[n]) + Im{y n (9)}sin(9[n]), (9) 

where Re{#} extracts the (real) I component of the symbol. Equation (9) is an estimation of the 
derivative d(y n (9)e' j9[n] )/d9 under the condition that y n (9)=y n . 

For the second configuration of FECL module 502a shown in FIG. 6, the data signal is 
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the output signal from combiner 602. The phase error term e S A-CM[n] of SA-CM detector 801 is 
computed with the data signal y n (6)e' je[n] -w[n]. The contribution D[n] is same as that in equation 
(9) given for the first configuration. 

For the third configuration of FECL module 502a shown in FIG. 6, the data signal which 
is the input of the e S A-CM[n] equation is the output decision d[n] from decision device 503. In the 
phase error term e SA -CM[n] of SA-CM detector 801 the variable y n (0) is replaced with the 
decision d[n] as given in equation (10): 



where f (•) is a non-linear function which stands for the decision function. The decision device 
as described in equation (10) operates on a real symbol. If we assume that f(«) is a differentiable 
function, then its derivative with respect to the variable "•" inside the brackets is denoted by 
f (•). The contribution DT[n] is then expressed as in equation (1 1), 



where Im{ y n (0)e" j9[n] -w[n]} in the expression of DT[n] is the same as in the first and second 
configurations. For a given implementation of DT[n], the value of the derivative of fi», shown 
as f (•), is approximated using the difference of adjacent symbols. 

As described previously, carrier recovery may also be implemented with the equalizer 
feedback filter operating in the passband. Returning to FIG. 7, FECL module 502b is shown 
including feedback filter 505 operating with input samples at the passband. The complex (FIR) 
signal (sample stream from forward filter 501) is applied to complex combiner 703. Complex 
combiner 703 combines I and Q components of the complex (FIR) signal with corresponding I 
and Q components generated by feedback filter 505. The input to feedback filter 505 is 
controlled by switch 704. When switch 704 is in position #1 the output of complex combiner 
703 is applied to feedback filter 505 (linear IIR equalizer operation). When switch 704 is in 
position #2, the output of complex combiner 703 is applied to phase de-rotator 701 and the 
output of phase re-rotator 702 is applied to feedback filter 505 (DFE operation). 

Phase de-rotator 701 and phase re-rotator 702 receive an estimate of 0[n] from carrier 
tracking loop 503, where 0[n] is an estimate at discrete time n of the phase offset. Phase de- 



d[n] = f(Re{y n (e)e- j9M -w[n]}), 



(10) 



DT[n] = f (Re{y n (e)e- je[n] -w[n]}) Im{ y n (0)e- je[n] -w[n]} 



(11) 
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rotator 701 operates on equalized complex (I and Q) components denoted y n (6)=I n (e) + jQ n (G), 
and computes the product y n (0)e- je[n] . Decision device 504 maps I and Q components of y n (9)e 
j9[n] to the nearest symbol d[n]. The symbol d[n] from decision device 504 is applied to phase re- 
rotater 702, which computes the product d[n]e +je[n] . Thus, for DFE operation, the feedback 
signal from feedback filter 505 is re-rotated to align in phase with the output signal from forward 
filter 501. 

Carrier tracking loop 503 generates an estimate of 0[n] provided to phase de-rotator 701 
and phase re-rotator 702. Embodiments of the present invention operating in the passband (e.g. 
FECL module 502b such as shown in FIG. 7) may generate an estimate of 9[n] using one of at 
least two configurations. These configurations provide an advantage in that both forward filter 
501 and feedback filter 505 process data that is not precisely at baseband. 

In the fourth exemplary configuration of FECL module 502b, carrier tracking loop 503 
receives Re{z[n]e" j6[n] }. The signal z[n] denotes the complex output signal from complex 
combiner 703 defined by z[n]= y„(0)-w[n], where y„(0) is the output of the forward filter 501 and 
w[n] is the output of the feedback filter 505. The value for e SA -CM[n] is computed with data 
signal z[n]e jew , where z[n]= y n (0)-w[n]. The contribution DT[n] for the fourth configuration is 
as given in equation (12): 

DT[n] = Re{y„(9)-w[n]}cos(0[n]) + Im{y n (9)-w[n]}sin(0[n]). (12) 

In the fifth exemplary configuration of FECL module 502b, the carrier tracking loop 503 
operates on the Re{d[n]}, where d[n] is the complex output of decision device 504. The value for 
e S A-CM[n] is thus computed with d[n]. The contribution DT[n] for the fifth configuration is as 
given in equation (13): 

DT[n] = f (Re{(y n (0)-w[nj) e - j9[n] }) lm{(y n (0)-w[n]) e - je[n] } (13) 
where the function f (•) is as defined above. 

FIG. 9 shows an exemplary implementation for portion 807 of carrier tracking loop 503 
operating in accordance with the first configuration of FECL module 502a of FIG. 6. At discrete 
time n, SA-CM phase detector 801 operates on the previous estimate of 0[n-l] and on the inputs 
Re{y n (0)} and lm{y n (0)}. The signal Re{y n [0]} is used to computed e SA . C M[n] by squaring 
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Re{y n (0)}, subtracting p 2 and multiplying the difference by Re{y n (0)}. The signal DT[n] is 
formed by i) multiplying the sine of 0[n-l] with Re{y n (0)}, ii) multiplying the cosine of G[n-1] 
with lm{y n (0)}, and iii) adding the two intermediary products of i) and ii). The phase detector 
signal is formed by multiplying esA-cMM with DT[n] to drive loop filter 802. Samples at 
5 discrete time (n-1) are employed in the calculation of DT[n] to specify that the current value of 
0[n] is derived from the previous estimate of 0[n-l]. 

Implementations for the second and third configurations of FECL module 502a, and for 
the fourth and fifth exemplary configurations of FECL module 502b, are derived in a similar 
manner to that shown and described for the first configuration of FIG. 9. One skilled in the art 
10 would realize that the implementation of FIG. 9 is exemplary only, and many different 
implementations may be designed that perform the calculations as outlined herein. 

FIG. 10 shows an implementation of module 807 when csa-cm[u] is calculated externally 
P to the carrier tracking loop 503 , such as when the receiver uses the SA-CM error term generated 
43 for adaptation of the taps of the equalizer. The value used for esA-CM[n] is the computed equalizer 
s _^5 error term csa-cmMeq (for example, from a look up table) which is also applied to SA-CM phase 
O detector 801. Otherwise, carrier recovery for the implementation of module 807 of FIG. 10 
S operates in a similar manner to that described with respect to FIG. 9. Thus, the processes of 
adaptively generating tap coefficients for equalization and estimation of carrier phase may be 
00 jointly optimized. 

&f) Some embodiments of the present invention may employ soft-switching criteria to 

U adaptively switch modes between SA-CM (blind) phase detection and a decision-directed (DD) 
phase detection for carrier recovery. Such adaptive switching may be based on relative quality of 
the received signal. As shown in FIG. 1 1, carrier tracking loop 503 may include signal quality 
measure (SQM) processor 1 103, phase detector switch 1 102, and DD phase detector 1 101 in 

25 addition to the SA-CM phase detector 801. DD phase detector 1101 determines a corresponding 
estimate of carrier phase error based on a cost-based error term. For example, DD phase detector 
1101 may employ a least mean square (LMS) error term ClmsM, when operating in a decision- 
directed mode. Other types of phase detectors may be further included that generate carrier 
phase error terms in accordance with other methods. SQM processor 1 103 may monitor the 

30 received signal quality to determine if estimates of the symbol value are relatively accurate, 
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based on some predetermined threshold. SQM processor 1101 may thus generate a signal 
quality measure (SQM) that is provided to phase detector switch 1 102. Using SQM, phase 
detector switch 1 102 selects one or more of the cost-based timing errors (e.g., esA-cM[n] or 
ClmsM). The selected cost-based timing errors may then be used alone or in combination by 
5 phase detector switch 1 102 for generating the estimate of angle 0. 

Signal quality measured by signal quality processor 1 103 may include: 1) signal-to-noise 
ratio (SNR) for the signal output from equalizer 461 (FIG. 4); and 2) SNR for the signal from, 
for example, a trellis decoder of detector 424 or decoder 425. In addition, signal quality 
processor 1 103 may use criteria for symbol-by-symbol switching between SA-CMA and LMS 

10 update modes for an adaptive equalizer. Such switching may be as described in a U.S. patent 
application serial no. 09/549,368, entitled "A HYBRID SOFT AND HARD DECISION 
FEEDBACK EQUALIZER," filed on April 14, 2000, which is incorporated herein by reference. 

^ Signal quality processor 614 may also use a prescribed number of symbol intervals between 

ffl signal (timing) acquisition and the beginning of symbol processing. Signal quality processor 

fjfo 1103 may use these factors alone or in combination to generate SQM. 

D One skilled in the art would understand that the equations described herein, such as 

J; equations (3) through (13), may include scaling, change of sign, or similar constant linear 

^ modifications that are not shown for simplicity. One skilled in the art would realize that such 

ofj modifications are readily determined or derived for the particular implementation. Thus, each of 

jSp the described equations may be subject to such modifications, and are not limited to the exact 

Q forms presented herein. 

To reduce stochastic jitter, error terms derived in accordance with the CM criterion may 
be modified by jointly equalizing and decoding. US patent application serial no. 9/099,730 
entitled "A METHOD OF ESTIMATING TRELLIS ENCODED SYMBOLS UTILIZING 

25 SIMPLIFIED TRELLIS DECODING," filed on November 20, 2000, and US patent application 
entitled "DATA SLICER FOR COMBINED TRELLIS DECODING AND EQUALIZATION," 
filed on November 20, 2000, describe various decision devices which perform partial trellis 
decoding and are incorporated herein by reference. For a VSB-modulated signal in accordance 
with the ATSC standard (an "ATSC 8-VSB signal"), the 8-level signal is decomposed into 

30 disjoint sets, each set with less than 8 elements (a "set-partitioned signal"). US patent 
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application serial no. 09/457,496, entitled "EQUALIZER METHOD AND APPARATUS 
USING BLIND CM A EQUALIZATION AND PARTIAL DECODING;' filed October 9, 1999, 
and incorporated herein by reference, employs the set-partitioned signals to define a new CM 
error term that incurs less stochastic jitter than the original CM error term. It is understood by 
5 one skilled in the art that such techniques to generate modified CM error terms may be employed 
in forming error signals in accordance with the present invention. 

Since the CM error term is a cubic function of the signal y n (9), particular 
implementations of the present invention using finite precision may require large bit- width 
multipliers and accumulators. To reduce the computational burden, low-complexity 
1 0 implementations of the CMA error term have been developed, such as those described in US 
patent application serial no. 09/231,888, entitled "ADAPTIVE EQUALIZER WITH 
ENHANCED ERROR QUANTIZATION," filed Jan. 14, 1999, and US patent application serial 
no. 09/261,843 "ADAPTIVE EQUALIZER WITH DECISION DIRECTED CONSTANT 
5 MODULUS ALGORITHM," filed March 1999, whose teachings are incorporated herein by 
35 reference. As would be apparent to one skilled in the art, the various methods of carrier recovery 
^ and synchronization described herein may be extended to include the use of quantized versions 
p of blind error terms. 

s Although the present invention is described herein with respect to the S A-CM cost 

rk function, the present invention is not so limited and may be extended in general to carrier 
8) recovery based on single-axis blind cost criteria as a class. For example, one skilled in the art 
O may extend the teachings herein to include criteria of Bussgang-class cost functions that may be 
related to the CM and SA-CM cost function described herein. The Bussgang-class of functions 
are well-known in the art and are described in, for example, Simon Haykin, Blind 
Deconvolution , Chapter 2, PTR Prentice-Hall, Englewood Cliffs, NJ (1994), whose teachings are 
25 incorporated herein by reference. 

The present invention has been described with respect to a digital television receiver, and 
in particular, to digital television signals transmitted as symbols with quadrature amplitude 
modulation (QAM), an phase-shift keyed modulation (PSK), or a vestigial sideband modulation 
(VSB). However, one skilled in the art would realize that the techniques described herein may 
30 be applied to any receiver processing QAM, PSK, or VSB signals. 
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As would be apparent to one skilled in the art, the various functions of carrier recovery, 
tracking, and carrier phase detection may be implemented with circuit elements or may also be 
implemented in the digital domain as processing steps in a software program. Such software 
may be employed in, for example, a digital signal processor, micro-controller, or general- 
purpose computer. 

The present invention can be embodied in the form of methods and apparatuses for 
practicing those methods. The present invention can also be embodied in the form of program 
code embodied in tangible media, such as floppy diskettes, CD-ROMs, hard drives, or any other 
machine-readable storage medium, wherein, when the program code is loaded into and executed 
by a machine, such as a computer, the machine becomes an apparatus for practicing the 
invention. The present invention can also be embodied in the form of program code, for 
example, whether stored in a storage medium, loaded into and/or executed by a machine, or 
transmitted over some transmission medium, such as over electrical wiring or cabling, through 
fiber optics, or via electromagnetic radiation, wherein, when the program code is loaded into and 
executed by a machine, such as a computer, the machine becomes an apparatus for practicing the 
invention. When implemented on a general-purpose processor, the program code segments 
combine with the processor to provide a unique device that operates analogously to specific logic 
circuits. 

It will be further understood that various changes in the details, materials, and 
arrangements of the parts which have been described and illustrated in order to explain the 
nature of this invention may be made by those skilled in the art without departing from the 
principle and scope of the invention as expressed in the following claims. 
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